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ignal propagation is essentially line of 

sight at the 500—3000-MHz frequencies 

being considered for satellite digital 
audio broadcasting (DAB). It is easy to design a 
satellite link when the receiver has a clear line of 
sight to the satellite, but very difficult when the 
view of the satellite is blocked by an obstacle such 
as a tree or building. 

The mobile reception case is in between 
these two extremes. Line-of-sight operation is inter- 
spersed by periods of signal blockage as the vehicle 
travels down the road. The times of blockage range 
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The Satellite Communications (SATCOM) Program is managed by the Jet Propulsion Laboratory for 
the National Aeronautics and Space Administration. The goals of this Program are to develop advanced 
system concepts and technologies in the areas of mobile, micro and personal communications. A major 
thrust of the Program through the mid-1990s is the development of the mobile terminal to operate with 
NASA’s Advanced Communications Technology Satellite. 


from a few milliseconds to many seconds. If the 
effects of at least the shorter signal outages can be 
eliminated by the receiver, satellite-only DAB re- 
ception may be practical over significant portions 
of the coverage area. 

Satellite DAB systems intended for urban 
and other difficult reception areas will generally 
require terrestrial booster transmitters to cover 
areas of difficult satellite reception. The experi- 
mental system discussed here does not include this 
feature; it is strictly a line-of-sight system. 


System Configuration 

The configuration of the Direct Broadcast 
Satellite-—Radio (DBS—R) experimental system is 
illustrated in Figure 1. The major components of 
the DBS—R experimental system are the modulator 
and demodulator, two Pacific Communication 
Sciences, Inc. (PCSI) codecs and the JPL Mobile 
Satellite Experiment (MSAT-X) transmitter and 
receiver. INMARSAT provided time on the 
MARECS-B2 Atlantic West satellite, a geostation- 
ary satellite located at 55° W. The limited power 
available from this satellite constrained the experi- 
mental data rate to below 20 kbps. 

The input audio, limited to a 3.5-kHz band- 
width, is digitized and compressed to a data rate of 
either 19.2 or 16.0 kbps by the PCSI codec, as 
chosen by the user. The codec puts out serial data 
in frames of 15 msec. Each frame contains a syn- 
chronization bit as well as some error correction 
coding. 

The DBS—R modem design, based on simu- 
lation results of [1] and previous JPL modem de- 
signs [2], is a fully digital modulator that supports 
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Figure 1. 
DBS-R system 
setup. 
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data rates of 19.2 and 16.0 kbps. The modulation 
utilized is digital phase-shift keyed with 100% 

root raised cosine pulse shape filtering for spectral 
efficiency. The coding employed is the optimal 

rate 1/2, K = 7, convolutional code with user select- 
able block interleaving of a full-size interleaver 

(16 x 2 K symbols), a half-size interleaver (16 x 1 K 
symbols) or no interleaver. 

The MSAT-X transmitter and receiver 
were previously used during the MSAT-X field 
experiments. Their main function is the up- and 
down-conversion to the desired frequencies of 
INMARSAT’s MARECS-B2 satellite (uplink 
C-band, downlink L-band). 

The DBS—R demodulator consists of three 
major components: an analog-to-digital (A/D) 
board, a digital front end board and a multiproces- 
sor demodulator board. The A/D board performs 
12-bit conversion at a sampling rate of 64 times the 
data rate. The digital front end board’s function is 
matched filtering. The multiprocessor demodulator 
board consists of five TMS320C25 digital signal 
processors and a Viterbi decoder chip. The majority 
of the demodulator processing is performed by this 
board (symbol timing algorithm, differential decod- 
ing, de-interleaving and 3-bit soft decision Viterbi 
decoding). In addition, the multiprocessor board 
provides diagnostic data (symbol error rate, power 


The symbol timing algorithm employed in the 
DBS-R modem proved to be highly robust ina 
fading environment. 


measurements and interleaver synchronization in- 
formation) that were recorded during the DBS-R 
field experiment. 


The PCSI codec decompresses the digital data 
from the modem and converts them back into an 
analog audio signal. The codec is self-synchroniz- 
ing but loses synchronization when a sufficient 
number of frame synchronization bits are in error 
— in which case it shuts down the audio output. 

If the data error rate suddenly becomes high, the 
codec can put out garbled audio before it enters the 
synchronization loss condition. When good data are 
restored, the codec resynchronizes and re-enables 
its audio output within approximately 250 msec. 

Differential encoding and differentially coher- 
ent detection were used in the DBS—R modem be- 
cause of the fading environment present in the 
mobile communications systems. Closed-loop sys- 
tems are more susceptible to outages due to fades 
and also take longer to recover from a fade than a 
typical open-loop system. The symbol timing algo- 
rithm employed in the DBS—R modem is shown in 
Figure 2. This scheme was also chosen because of 
earlier successes during the MSAT-X experiments, 
where it proved to be highly robust in a fading en- 
vironment. The random-walk filter in this process 
acts as a nonlinear lowpass filter that smooths out 
the symbol timing jitter. The chosen value of the 
random-walk filter should be large enough so that it 
can withstand moderate fades, but not so large as to 
hinder the symbol timing algorithm’s ability to up- 
date the receive clock due to any differences in the 
two clocks. 

The generator polynomial of the particular 
rate 1/2, constraint length 7, convolutional code is 
171,133 (expressed in octal). This code is the opti- 
mal d,.. (d,,,.= 10) code and is a standard commer- 
cial and scientific code used for space and satellite 
communications. Thus, there are many single-chip 
implementations commercially available, allowing 
for some degree of freedom in the design of the mo- 
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dem. Any convolutional code of constraint length 
greater than 7 was not considered feasible for the 
operational data rates of the modem. 

The function of the interleaver is to decorre- 
late error bursts to increase the Viterbi decoder’s 
ability to correct the errors. A 1-sec limitation on 
the size of the interleaver was deemed necessary for 
rapid reacquisition of interleaver synchronization in 
a fading environment. For this reason, an overall 
size of 32 K symbols was chosen as the maximum 
interleaver size (approximately 0.85 sec). The 
16 x 2 K symbol dimension was chosen because of 
the ease of modem software implementation since 
the TMS320C25 digital signal processors have a 
16-bit word length. 

In addition to the selection of the interleaver 
dimensions, an interleaver synchronization algo- 
rithm had to be developed. Each interleaver had 
two synchronization words located at the positions 
shown in Figure 3. These two words are each a 
127-symbol PN sequence and its complement. 
These two unique words were chosen because of 
the properties of PN sequences to correlate well 
with shifted versions of themselves. These two 
unique words are placed in the interleaver by writ- 
ing over the data symbols. This creates an 8/15 
punctured convolutional code for a small percent- 
age (12.4%) of the data. Simulation results have 
shown the degradation between the rate 1/2 and the 
punctured rate 8/15 convolutional codes to be negli- 
gible [approximately 0.1 dB over the bit error rates 
(BERs) of interest]. 

In order for the demodulator to declare 
interleaver synchronization, two consecutive syn- 
chronization words must be detected. The detection 
mechanism used in the demodulator is a bit-by-bit 
correlation between the received data signal and a 
known copy of the PN sequence. When a threshold 
level of 100 or more like-bit comparisons (or fewer 
than 27 like-bit comparisons for the complement 
synchronization word) is reached, detection of a 
unique word is declared. The choice of two con- 
secutive synchronization words of length 127 sym- 
bols — as opposed to one synchronization word of 
length 255 symbols — was made for two reasons. 
First, the processing complexity of detecting a syn- 
chronization word of length 127 symbols is half as 
much as for a word of length 255 symbols. This 
was a major consideration for the modem’s opera- 
tional data rates. Second, because of the fading 
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DBS-R symbol 
timing and 
differential 
decoding 
algorithm. 
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environment, the possibility of falsely declaring 
synchronization is more likely in the single syn- 
chronization word algorithm than in the chosen 
algorithm. Another reason for choosing four 
consecutive missed synchronization words is to 
provide some temporal diversity for the synchroni- 
zation algorithm in the fading environment. 

The demodulator checks for data integrity by 
monitoring two synchronization words embedded 
in each data interleave. When the demodulator has 
declared interleaver synchronization, data are trans- 
mitted from the demodulator to the PCSI codec. 
When four consecutive synchronization words go 
undetected, the demodulator is out of interleaver 
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Figure 3. 
Synchronization 
word locations 
in interleaver. 
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synchronization and discontinues data flow to the 
codec until resynchronization occurs. This feature 
can also be used to shut down the system during 
fades that are longer than two interleaves. Experi- 
mental results have shown, however, that the codec 
tended to lose synchronization before an out-of- 
sync condition was declared by the interleaver. 
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Figure 4. 
Fading channel 
simulator test 
configuration. 


dnterleaver Performance 

When a receiver is mobile, the link suffers 
from periodic shadowing from trees, buildings and 
other obstructions. Signal blockages generally re- 
sult in very deep fades, which essentially result in a 
BER of 0.5 during the period of the fade. [If the 


interleaver span and depth are Nand N’ dep SYD 


bols, respectively, and the fade duration is NV... sym- 


fade 
bols, the average error rate over the full interleaver 


is Nad (N pene Na or approximately 0.5.] 

One of three conditions may result, depend- 
ing on the length of the fade. If the fade is short 
enough, the interleaver will keep the average error 
rate to the Viterbi decoder low enough for it to cor- 
rect most errors. If the fade is longer, the average 
error rate will become too high for the decoder to 
correct, and eventually the voice codec will lose 
synchronization. If the fade lasts for longer than 
two interleaves, the modem will declare interleaver 
synchronization loss and will turn off the data to the 
voice codec. 


Laboratory Test Results 

Prior to the DBS-R field experiment, the 
modem was subjected to laboratory tests which 
simulated a Gaussian channel and three different 
types of fading channels. The purpose of these tests 
was to define the approximate modem performance 
for the field experiment. The laboratory test was 
set up to run at a carrier frequency of approximately 
50 MHz with additive white Gaussian noise 
(AWGN) and fade effects, as shown in Figure 4. 
The fade data used in these simulations were ob- 
tained during a previous field experiment; the sta- 
tistical information was obtained from [3]. The 
complete fade distribution for these three data files 
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is shown in Figure 5(a). The fades have been cate- 
gorized as light, moderate or heavy. The probability 
of a fade greater than 3 dB for the light fading case 
was about 8%. For the moderate fading environ- 
ment, the probability of a fade greater than 3 dB 
was 12%, and for the heavy fading channel this 
value was 45%. BER results were obtained for all 
three of these fading environments for E,/N_ rang- 
ing from 7 to 12 dB. These results are shown in 
Figures 5(b)}-(d). At a BER of 10°, the gain of the 
full interleaver over no interleaver for the light fad- 
ing environment was approximately 0.8 dB. For the 
moderate and heavy fading environments, the gains 
were 1.3 and 1.6 dB, respectively. 


Field Experiment and Laboratory Test Parameters 

The focus of the DBS-R field experiment 
was the performance of the different size interleav- 
ers in various mobile environments while varying 
the available link margins. The experiment was 
conducted on different segments of highway in 
Connecticut. These stretches of road were chosen 
on the basis of varying degrees of satellite signal 
blockage. For each of these stretches of road, three 
identical experimental runs were made: one using 
the full-size interleaver, one using the half-size 
interleaver and one with no interleaver. The trans- 
mitted signal power was also varied by several 
decibels to test the performance of the system at 
different signal-to-noise ratios. The BERs for each 
of these setups were recorded and will be compared 
to each other in the quantitative results section of 
this article to determine the performance gain as the 
interleaver size is increased. The important param- 
eters in this comparison are the interleaver size, 
BER, E,/N_, fade depth and duration of fade. A sec- 
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ondary comparison to be made involves how well 
the interleaver synchronization algorithm performs 
for each interleaver size. The key parameters for 
this setup are interleaver synchronization losses, 
E,/N_, fade depth and duration of fade. 


Experiment Setup 

The setup for the experiment is shown in Fig- 
ure 6. The broadcasted signal was simultaneously 
received at COMSAT’s Earth station and in the 
DBS-R experimental van. The main difference be- 
tween the two setups can be seen in the antennas. 
The Earth station had a high gain stationary antenna 
with a gain of about 24.0 dB, and the experimental 
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van had a medium gain, tracking (Yagi) antenna 
with a gain of 12.25 dB. Since the gain of the sta- 
tionary antenna at the Earth station was so much 
higher than the Yagi antenna, the signal captured 
was practically error free (the system BER was less 
than 10°). Since the satellite signal power tended to 
vary with time, the data recorded in the Earth sta- 
tion were used as a reference for the signal strength 
data recorded in the DBS—R experimental van. 

The data recorded during the experiment at 
the Earth station were pilot signal level, BER infor- 
mation, unique word detection count, interleaver 
synchronization loss count, data signal power level, 
time and audio signal. In the DBS—R experimental 
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Figure 6. 
DBS-R field 
experiment 
test setup. 
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van, the data recorded during the experiment were 
pilot signal level, BER information, unique word 
detection count interleaver synchronization loss 
count, data signal power level, Loran-C information 
regarding the vehicle location, antenna pointing in- 
formation (azimuth and antenna tracking informa- 
tion), time, audio signal and a video recording of 
the line of sight between the van and the satellite. 
Most of this information in both the COMSAT 
Earth station and the DBS-R experimental van was 
recorded using the MSAT-X Data Acquisition Sys- 
tem, a dedicated personal computer interfaced to 
two Bernoulli cartridge drives for the sole purpose 
of mass storage of data. 


Field Experiment Link Budget 

Table 1 outlines the best case link budget for 
the DBS-R field experiment for the two operational 
data rates of the modem. The link budget shown in- 
volves the COMSAT Earth station as the transmit 
side of the link; INMARSAT’s MARECS-B2 sat- 
ellite, and the DBS—R experimental van as the re- 
ceive side of the link. The actual satellite power 
available (determined during previous experiments 
[4]) over the link was at times up to 3 dB lower 
than stated in Table 1, so that the link margin 
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went down to about 1.7 and 2.5 dB for 19.2 and 
16.0 kbps, respectively. The fluctuation in the avail- 
able satellite power occurs because the satellite 
transponder is equipped with an automatic level 
control to maintain the total satellite output at a con- 
stant level of 34.50 dBW. Also, it should be noted 
that there was a varying amount of other traffic on 
the satellite. 


DBS-R Field Experiment 

The DBS-R field experiment was performed 
from September 7 to September 24, 1991. On Sep- 
tember 7, the DBS—R experimental van departed 
Los Angeles for Southbury, Connecticut. For 
3 hours each day (from 3 p.m. to 6 p.m. Eastern 
Daylight Time), the DBS—R experimental system 
was operational. The reason for running the experi- 
ment while the van was en route to Connecticut was 
to observe the effects of different elevation angles 
of the MARECS-B2 satellite. The Yagi antenna 
was able to track the DBS-R pilot signal at an ele- 
vation angle as small as 12 deg (Los Angeles). 
However, the antenna gain did not get high enough 
to be able to consistently record audio until the 
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Parameter Modem data Modem data 
rate = 19.2 kbps rate = 16.0 kbps 

Available power, dBW 34.50 34.50 
% Satellite power, dBW 18.18 18.18 
Transmitter 

Power, W 4.00 4.00 

Power, dBW 6.02 6.02 

Antenna gain, dB PAU 21.18 

EIRP, dBW 25.00 25.00 
Path 

Space loss, dB —187.80 —187.80 
(se lel ll bs eae Pace esrtnee sleat teste aan te ie a ate 
Receiver 

Antenna gain, dB 12.25 1225 

Pointing loss, dB 0.00 0.00 

Signal power, dB —153.60 —153.60 

System temperature, dBK 24.50 24.50 

N_., dBW/Hz —204.10 —204.10 
Channel 

Received C/N ,, dB 5550 33.95 

Link E,/N_,, dB 10.72 | Rew 

E,/N, required, dB 6.25 6.25 

Implementation loss, dB 1.00 1.00 
8 Das ee Re 
Available margin, dB 4.47 5.26 
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Table 1. 
L-band DBS-R 
experiment link 
budget. 
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Figure 7. 
Route of 
DBS-R 
experimental 
van from 

Los Angeles, 
California, 

to Southbury, 
Connecticut. 


Table 2. 

Yagi antenna 
gain for various 
elevation 
angles. 
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elevation angle was at least 22 deg (just west of 
Oklahoma City, Oklahoma) [5]. In traveling from 
Los Angeles to Southbury, the antenna elevation 
angle varied from 12 to 39 deg. Figure 7 shows the 
DBS-R experimental van’s route across the United 
States and Table 2 shows how the Yagi antenna’s 
gain varied with different elevation angles. 

While the van was traveling to Connecticut, 
the only mode in which the DBS—R modem oper- 
ated was a data rate of 16.0 kbps with a full-size 
interleaver. This allowed the system every possible 
gain to compensate for the lower look angles. Once 
the van reached Connecticut, all modes of operation 
were tested. The main objective of the experiment 
— to test the performance of the DBS—R modem 


Elevation angle, deg Antenna gain, dB 


5 7.00 
10 furs 
15 9:25 
20 10.00 
25 10.25 
30 1125 
Ey) 12.00 
40 25 


Experiment Held in 
Southbury, Connecticut; 
Elevation Angle = 39 deg 


Start of Coherent Audio for 
DBS-R Field Experiment; 
Elevation Angle = 22 deg 


with different interleaver sizes in various fading en- 
vironments — was accomplished in Connecticut. 
Several different stretches of roads that exemplified 
different types of fading channels were chosen. 
Three different runs were made on each of these 
roads — one for each different size of interleaver. 
These tests were performed from September 16 to 
September 22. 


Qualitative Results 

When the DBS-R system in the van was ina 
clear line of sight, there was no difference between 
the audio quality in the van and in the COMSAT 
Earth station. Due to the bandlimiting of the signal 
(3.5 kHz), the signal did not have the quality of an 
FM broadcast, but the signal quality was compa- 
rable to an AM broadcast. The full-size interleaver 
helped the signal quality in the fading environment, 
as many of the shorter fades that affected the no- 
interleaver configuration were eliminated. How- 
ever, even in the lightest fading test runs, there were 
some stretches of road where the signal blockage 
events were too long for even the full-size inter- 
leaver to combat. When outages did occur, the 
PCSI codec lost synchronization but occasionally 
put out a short burst of garbled audio before becom- 
ing silent. A possible solution to this problem is 
presented below. 

One of the drawbacks of the experimental 
system was that the voice codec ran independently 
from the demodulator. Thus, voice codec fades 
were not controlled as well as they could have been. 
The Viterbi decoder chip has an uncoded symbol 
error estimation feature which can be used to moni- 
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tor the quality of the link [6]. The uncoded symbol 
error estimator operates by re-encoding the Viterbi- 
decoded data and then comparing this data sequence 
to the received, predecoded data sequence. The 
number of errors in a user-specified block size is 
then accumulated. This symbol error estimator was 
very useful in determining system performance dur- 
ing the field experiment and is a feature that will be 
used for codec muting control in further work on 
this system (see Figure 8). 

As seen in experimental testing, the codec 
loses synchronization when the error rate is above 
a certain level (uncoded symbol error rate = 
6.0 x 10°’, corresponding to a coded bit error rate 
of 2.0 x 10°). Below this threshold, the audio qual- 
ity was acceptable and did not suffer gradual degra- 
dation. Using the symbol error estimator to mute 
the audio during periods of high errors will result 
in a higher quality audio output (see Figure 9). 


Quantitative Results 

Figure 10(a) categorizes each of the three 
fading environments in terms of their fade distribu- 
tions. As in the laboratory simulation results, these 
three fading environments are defined as light, 
moderate or heavy. The probability of a fade greater 
than 3 dB for the light fading environment was 
9.5%. For the moderate fading environment, the 
probability of a fade greater than 3 dB was 15% and 
for the heavy fading channel this value was 37%. 
Figures 10(b)—(d) provide the corresponding BERs 
for each of these fading environments for the full- 
size interleaver, half-size interleaver and no-inter- 
leaver modes of operation. Table 3 summarizes the 
BER results for the different fading environments 
for the full-size interleaver mode of operation. The 
full-size interleaver case provides noticeable im- 
provement over the no-interleaver case, as outlined 
in Table 4. Corresponding results for the half-size 
interleaver case are shown in Tables 5 and 6. 

Analyzing the data pertaining to the interleav- 
er synchronization algorithm showed that the mo- 
dem maintained synchronization in fades that were 
less than about 2.2 sec for the full-size interleaver 
case and in fades that were less than about 1.0 sec 
for the half-size interleaver case. These results were 
about as expected — in order to lose interleaver 
synchronization, two complete interleaver’s worth 
of data must contain fade data. This corresponds to 
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withstand can be explained by the fact that if the 
fade does not occur exactly at the beginning of the 
interleaver, the algorithm becomes even more ro- 
bust in fades. 
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Conclusion 

The DBS-R field experiment has shown that 
the use of an interleaver significantly improves the 
performance of the DBS—R modem in a fading en- 
vironment (by as much as 1.3 dB in certain fading 
environments). However, coding and the use of the 
interleaver will not fully combat all the problems 
that are encountered in a mobile environment. 
There were still significant signal outages even in 
the most favorable of the fading environments. 

An addition to the link margin will help to 
improve system performance. However, further re- 
search and development are necessary to investigate 
other fade mitigation techniques. 

At the present time, further development on 
Direct Broadcast Satellite Audio is under way at the 
Jet Propulsion Laboratory. The main objective of 
this work is to devise new techniques to improve 
system performance in difficult reception environ- 
ments. The techniques being investigated through 
analysis and simulation include further refinement 
of the interleaver size and dimensions, as well as 
channel equalization and spatial diversity techniques. 
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hen the commercial use of satellites 

for telecommunications began, 

C-band (6/4 GHz) was the allocated 
frequency range. The limited bandwidth (500 MHz) 
allocated in C-band soon became congested. In 
the 1980s, a higher frequency band, Ku-band 
(14/12 GHz), was allocated specifically for do- 
mestic SATCOM use. The frequency allocation 
(500 MHz) in this band is also becoming saturated. 

This has led to the investigation of the next 
higher frequency band, Ka-band (30/20 GHz), to 
determine its utility for satellite communications. 
The wider bandwidth (in excess of 1000 MHz) 
available in this band is likely to be heavily used. 
However, the wavelengths in this band are more 
susceptible to perturbing elements of space-to- 
Earth communication links. These perturbing ele- 
ments are shown in Table 1. 

Rain attenuation is dependent upon rain rate, 
rain-path length and drop-size distributions. At 
Ka-band (1—1.5-cm wavelength), the drop size is 
comparable to the wavelength. Thus, attenuation is 
severe, even to the extent that in some cases it be- 
comes economically unattractive to provide com- 
pensation for rain fade. As a result, communication 
link outage is inevitable. 


For cost-effective business applications, the 
communications system should be physically com- 
pact. Antenna size should be limited to about 1 m 
or smaller in diameter. This small antenna size, 
coupled with the limited RF transmit power avail- 
ability from miniterminals, imposes a restriction 
on the range of the power control in the 30-GHz 
uplink. On the other hand, the 20-GHz downlink 
power from the satellite transponder is limited, and 
adaptive control of downlink power is difficult to 
manage in large networks covering different cli- 
matic regions. Normally, downlink fade is compen- 
sated for by a finite margin allocation. In some 
applications, such as the Advanced Communica- 
tions Technology Satellite (ACTS), the downlink 
margin is augmented by reducing the data rate and 
adaptive application of forward error correction 
techniques. For example, a data rate reduction by a 
factor of four enhances the power margin by 6 dB. 
Another 4 dB of margin can be introduced by rate 
(1/2) coding. These enhancements to the margin, in 
addition to the built-in downlink margin, will in 
most cases circumvent the 20-GHz downlink fade. 

Normally, the extent of the rain fade is mea- 
sured by comparing the downlink beacon signal to 
the clear-sky measurement. However, the clear-sky 
measurement is corrupted by gaseous attenuation, 
clouds, antenna pointing error, satellite movement 
and scintillation. These perturbations are measur- 
able and can be eliminated from the rain-fade effect 
by combining radiometric observations with bea- 
con measurements. The measured attenuation, 

A, (dB) at frequency f,, can be translated to atten- 
uation A, at f, by scaling law as follows: 


2 
A, (dB) = ae (dB) 


1 


For example, if A, = 3 dB at 20 GHz, then at 
30 GHz 
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In some applications (ACTS, for example), 
beacons are provided at both uplink and downlink 
frequency bands, so independent signal fade mea- 
surements can be performed. 

The basic system concept of rain-fade detec- 
tion is illustrated in Figure 1. The process control 
sequences for uplink rain fade compensation are il- 
lustrated in Figures 2(a) and 2(b). The system time 
lag comprises delays in filtering, measurement, sig- 
nal processing and control actuation. 

In essence, the problem is: Given a measured 
or calculated attenuation at time ¢, what is the ex- 
pected attenuation at (t +t), where t is the system 


| | Satellite 


/ 
3 Beacon at f, or 


Beacon 
Receiver 


f, and/or f, 


Beacons at f and f, 


time lag? 
Element Effect Remedy 
Gaseous attenuation Noticeable Built-in margin 
Hydrometer attenuation 
Cloud Measurable Built-in margin 
Rain Severe Adaptive compensation 
Snow Measurable Can be circumvented 
Fog Negligible N/A 
Depolarization Complex Can be circumvented 
Scintillation Rapid fluctuation of signal; Not easy to circumvent 


dependent on site, time 
of day and season 


Excess noise emission Signal absorption gives rise 
to broadband noise emission additional system margin 


Transmit 
Measurements Calculation Process Power 
at I, for f, (Uplink Power) 
by eee uh wate System Time Lag fs ea 


Control 
Measurements Signal Process 
at 9 (Uplink Power) 
—— System Time Lag omits che ean! 


Can be circumvented by 


Downlink 
Beacon at f, 


Control 
Signal 


Uplink 


Transmit 
Beacon at f, 


Power 


Figure 1. 
Basic system 
concept. 


Table 1. 
Perturbing 
elements of 
Ka-band space- 
to-ground com- 
munications. 


Figure 2(a). 
Case when 
only downlink 
beacon is used 
to predict 
uplink fade. 


Figure 2(b). 
Case when 
uplink beacon 
is available. 
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Figure 3. 
Uplink 

power control 
fundamentals. 


Mathematical Review 

We now discuss two basic axioms [1] of con- 
trol engineering for uplink power control as illus- 
trated in Figure 3, where Z(f) is the input signal 
sample, e.g., rain attenuation as a function of time; 
(t,t) is the impulse function of the signal transfor- 
mation process, and X(f) is the transformed control 
signal. It can be shown that 


t 


X(t) = if h(t-t)Z(t) dt (1) 
Input Z (t) Transformation Control X (t) 
h (t, t) 


Axiom 1 


The output of a matched filter is the convolu- 
tional integral 


u(t) = il h(t-t)r(t) dt (2) 


0 
at t= T (clock rate, for example). 


tip 
u(t) = if A(T-1)r(t) dt (3) 


0 


where /(¢) = filter impulse response and r(f) = 
receive waveform. 

The term u(t) is equivalent to the integrated 
output of the cross-correlation if 


A(T-t) = S(t) (4) 


where S(t) = correlation signal 


T 


ut) = f S(s)r(s)dt (5) 


0 


Equations (2-4) are important in the 
sense that we shall use these properties of cross- 
correlations in hardware design for uplink fade 
compensation. 


Axiom 2 


The rain rate, and the consequent attenuation, 
can be approximated by a random function of time. 
Therefore, exact values for the output signal, which 
result from input signal transformation at any mo- 
ment of time, cannot be given. Instead, the problem 
of specifying the output signal reduces to the prob- 
lem of specifying a system of functions that deter- 
mines the probability of the limits within which the 
value of the random signal exists at one moment of 
time and for a sequence of moments of time. 

Of particular significance in the theory of 
control engineering is a special class of random 
function — the class of stationary random functions 
(processes) wherein the time axis of arbitrary 
length t, for every moment of observation ¢, (k= 1, 
2, ..., 2), does not change the character of the distri- 
bution of the process. However, rain events cannot 
be strictly classified as a stationary random process. 
Perhaps it will be more appropriate to classify a 
rain event as a conditionally stationary random pro- 
cess. Under this condition, we can assume rain at- 
tenuation data samples derived from a rain event as 
a Markov process in which a sample is influenced 
only by the immediately preceding sample (first- 
order Markov process). 

One can find the properties of transforma- 
tions of random signals from their moment charac- 
teristics. Of special interest are the moments of first 
and second order, denoted by m_(t) and B (t, and 1,). 

The second moment of the unbiased random 
function of the process is 


R{t,,t) iz Bt, .t) ie m_(t,) m,(t,) (6) 


For stationary processes, ™m_(t) = m,constant, 
hence 


Rt,,t) es Rt,-t,) e Rt) (7) 


AP Rita 1992 


where t = (7, — 1). 

Equation (7) leads to the proposition that the 
correlation function for the random process Z(t) de- 
pends only on the difference of the arguments ¢, 
and t,. We shall use this important property of cor- 
relation function for the actual design of the uplink 
fade control in subsequent discussions. 


Input Signal Characteristics 

Under NASA funding and JPL sponsorship, 
Virginia Polytechnic Institute is carrying out radio- 
wave propagation experiments at 12.5, 20 and 
30 GHz using the European Space Agency’s Olym- 
pus Satellite, which is visible from Blacksburg, 
Virginia, at about 14-deg elevation. The first-year 
data collection/analysis report [2] has been pub- 
lished recently. This report analyzes the first year 
of rain attenuation and radiometric data. 

Rain attenuation is normally presented as a 
time series, a sample of which is shown in Figure 4. 
This figure also shows the contributions from gas- 
eous attenuation, diurnal fluctuations and scintilla- 
tions [2]. Simultaneous radiometric measurements 
are also made to measure sky brightness tempera- 
ture as shown in Figure 5 [2]. The sky brightness 
temperature increase due to signal absorption is 
governed by the relationship 


I, | (a-1) 
jee atten Mohs 
sky a a m 
therefore 
noe (8) 
leas 


where & = attenuation ratio due to rain absorption, 
T= cosmic background temperature = 2 K and 
T,, =rain medium temperature = 275-280 K. 
Therefore, one could accurately derive rain 
attenuation within a given range from sky bright- 
ness temperature measurements and then cross- 
check the beacon attenuation measurements exem- 
plified in Figure 4. When properly calibrated, the 
radiometer can measure gaseous attenuation in the 
clear-sky condition in an absolute sense. During 


Relative Power, dB 


0 3 6 ) 12 
GMT, Hours, October 22, 1990 


Sky Brightness Temperature, K 


0 
0.0 2.4 4.8 7.2 9.6 12.0 14.4 16.8 19.2 21.6 24.0 
GMT, Hours, October 22, 1990 


rain, radiometrically derived attenuation contains 
the contribution due to both gaseous and rain at- 
tenuation. Therefore, if the predicted attenuation 
exceeds a threshold set for the clear-sky condition, 
the hydrometer effects (attenuation due to rain) and 
the gaseous attenuation can be separated. On the 
other hand, beacon signals have a problem in set- 
ting an absolute reference level as well as problems 
with diurnal fluctuations. By proper calibration and 
manipulation, the effects due to gaseous attenua- 
tion, diurnal fluctuations and scintillation can be 
removed from Figure 4 to derive the net effect due 
to rain. 


Input Signal Dynamics 

The dynamics of rain-induced fade on satel- 
lite radio links is a complex subject that has been 
studied by various researchers [3-6]. The fade 
slope parameter (dB/sec) can be used to predict at- 
tenuation at time (t + t) when the sample value at 
time ¢ is known. 

Consider a radio path with a high gain an- 
tenna at each end. The majority of the energy trans- 
ferred between the two antennas passes through the 
volume defined by the first Fresnel zone. The first 


Figure 4. 
Beacon 
attenuation at 
12.5, 20 and 

30 GHz over a 
24-hour period. 


Figure 5. 

Sky noise 
temperature 

at 12.5, 20 and 
30 GHz over a 
24-hour period. 
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aE 


Figure 6. 
The first 
Fresnel zone. 


Fresnel zone is a prolate ellipsoid of revolution and 
is shown in Figure 6. The path attenuation is related 
to how much water is in the volume, but the rate at 
which the attenuation changes is related to how fast 
this volume fills. 


Path Length 


Antenna 


Projected Rain Area 


X 


An approximate mathematical relationship of 
fade slope in a terrestrial path in terms of rain rate 
and drop velocity [6] can be written as 


a 


ee 
G0) = 1 5{apR} v2) 1-4 (9) 
dt a 


where o = rain attenuation; a(f) and b(/) are 

the rain attenuation parameters; R = rain rate; 

v,= weighted average of drop velocity influenced 
by storm, wind velocity, etc.; b = major axis of 
the ellipsoid; a = minor axis of the ellipsoid, and 
t= time. 


for the fade slope occurs just at the beginning of the 
fade and when the attenuation reaches its maxi- 
mum. Thus, fade slope is not proportional to at- 
tenuation [6]. Furthermore, the rain drop velocity is 
influenced by a number of factors, including storms 
and wind velocity, and hence fade slope is not con- 
sidered a good measure for predicting uplink power 
control [6]. Fade as a function of time is perhaps a 
better measure for uplink power control. 


Control Algorithm Model 

Following Axioms 1 and 2, we propose the 
output signal manipulation. Given a set of se- 
quences, say d,, d,, ..., a, at ume t; Lyte what is 
the estimate of a at time t;? This set of sequences 
can be compared with attenuation sample values, 
e.g., a, (dB), a, (dB), ..., 4; (dB) at time ¢, (sec), 

t, (sec), ..., t,(sec), where (t,—¢,) = 1/m, t = system 
time lag and m is large but finite. The problem un- 
der consideration, then, is: If we can generate a pat- 
tern (history) of Den a,, how do we estimate 

a, (dB) at time t (sec)? 

This branch of science has found extensive 
applications in telecommunications engineering in 
the form of the Maximum Likelihood Sequence Es- 
timation (MLSE) [8]. The essence of MLSE is dis- 
cussed below as a prelude to the conception of an 
estimator model and its possible application in up- 
link fade prediction. 

Using a data communications analog, let us 
assume that the sequence of information amplitudes 
to be transmitted is donated by d,, d,, ..., 4, where 
nis large but finite. The purpose of MLSE is to op- 
erate on the receive signal R(t),0 <t<T"', where 


Maximum Likelihood Sequence Estimation (MLSE) can be used to develop a conception of an 
estimator model and its possible application in uplink fade prediction. 


Examination of Equation (9) reveals that the 
fade slope is nearly proportional to the rain rate R 
since b( f) = 1, but the slope varies with the third 
power of the drop velocity. The drop velocity de- 
pends on the drop size; hence, the assumption of v , 
as the weighted average of drop velocity due to dif- 
ferent drop sizes makes Equation (9) a more realis- 
tic solution. Drop velocities range from 2 to 9 m/sec 
for drop sizes of 0.05 to 0.7 cm [7]. Due to the na- 
ture of the Fresnel zone filling, the minimum value 


T! =nT + (s/m) — B = 1/T = channel bandwidth — 
to produce an estimate (4, d,, ..., @,) of the se- 
quence of transmitted signals (a,, d,, ..., @,), given 
that R(t) is perturbed by various transmission im- 
pairments. In the rain attenuation case, we have a 
similar situation, and d,, d,, ..., @, could be an actual 
rain attenuation data sequence, but the error in mea- 
surements and scaling law produced an estimate of 


ee ga Mego 3 Since the exact value cannot be mea- 


1? 
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SE ae 


sured with certainty, we attempt to minimize the 
probability of sequence error, i.e., the probability 
G,, G,, ..., Z, is different from d,, @,, ..., a, MLSE 
results in the minimum error probability. 

To define MLSE, the conditional probability 


density function must be defined first: 
P[R(t),0stsT",/a, =), Gg0 = 145.0, 80) 


T! 
K exp { - [ [R(t)-S()Pdt} (10) 
0 


0 


where R(t) = distorted receive signal, 


se) =f |S. ame-i0) 


i=l 


(the functional representation of the response of 
R(®) in the baseband form), K = constant and N, = 
noise density or system disturbance. 

The baseband receive signal S(t) due to input 
signal R(t) is analogous to X(t) due to Z(t) as dis- 
cussed above. 

The quantity J [+]? in Equation (10) is the 
mean square error (MSE), which is to be mini- 
mized. When the quadratic term in Equation (10) 
is expanded, we obtain 


T) T} T! 


if [R(t)]?dt-2 i R(t)S(t)dt + i} [S(t)]?dt 


0 0 0 


Finally, the objective function to be mini- 
mized is 


T! 
D[] = -2 if R(f)S(t)dt (11) 
0 


which is the correlation function discussed in Axi- 
oms | and 2. 

As discussed above, the correlation function 
for the conditionally stationary random process Z(t) 


depends only on the difference of the arguments ¢, 
and t,, or t, and (t, + t/m). Here, m is arbitrarily 
chosen for nearly a continuous check on the mea- 
sured parameter, e.g., rain attenuation as a function 
of time. 

Based upon these discussions, we can gener- 
ate the concept of an output signal estimator model 
for uplink fade control as shown in Figure 7, where 
the raw input beacon data are prefiltered to elimi- 
nate or manipulate scintillation, gaseous attenuation 
and diurnal effects. The filtered data are delayed 
by t seconds (the system time lag), and the de- 
layed signal is sampled with the sampling clock 
(rate = m) where m is to be decided by trial and er- 
ror. These sample data represent the sequence a,, 
Cried, AS discussed before. Another branch of 
prefiltered data passed through the system memory 
module, which modifies the input to account for the 
probable system perturbations resulting in the esti- 
mated sequence G,, a, ..., @, at the output of the 
sampler. The threshold level in the adaptor is set 
to bias out the effects of small-intensity scintilla- 
tions, gaseous attenuation and diurnal events. The 
adaptor output (4, G,, ..., d,) is now compared to the 
sampled sequence of the raw input d,, d., ..., @; in 
the MSE correlator. The MSE correlator output is 
further processed to generate the best estimate of 
the control signal. 


Practical Form of Rain Fade Compensation Estimator 

Based upon a first-order Markov process, a 
practical form of rain fade compensation estimator 
is proposed. 

A practical form of MLSE is an Adaptive 
Threshold Detection with Estimated Sequence 
(ATDES) [9]. In order to understand the basic prin- 
ciple of ATDES, we have to establish the principle 
of decision boundary shift as shown in Figure 8, 
where it is assumed that three consecutive samples 
of measurements are under examination: @_,, 4, and 
G,, at time (t—1/m), t and (t + t/m) respectively, 
where d_, represents the leading sample, d, is the 
objective sample and 4, is the trailing sample. Since 
the rain rate and the resulting attenuation changes 
with time, the decision boundary (threshold) due to 
G@_,,4, and 4, will differ and produce control errors 
from sample-to-sample observations. The problem 
is to compensate for the shift in the decision thresh- 
old from sample to sample. A practical scheme for 
this compensation is illustrated in Figure 9. 
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Decision Boundary (Threshold) 
for Actual Signal 
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+ 


Receive Samples Value 


Let the three consecutive receive samples be 
represented by 


by = flG_1,45,4,) + Ny (12) 


where N, is system disturbances. The optimum 
decision on J, is obtained by setting the decision 
threshold at d,, where 


evar A ; : 
dy 7 4 1,4,) a ECGwy-¥,a;)| (13) 


with a+1-V decision range, for example. 

To obtain the optimum decision threshold, 
ATDES first estimates the actual measured samples 
G_,, 4, and 4, by a tentative decision circuit operat- 
ing on the sign of input succeeded by a shift regis- 
ter. There, the suboptimum threshold d, for b, is 
obtained by 


d, a [F(@_,,1,4,) + F(@_,,-1,4,)] (14) 


Figure 9. 
Simplified 
block diagram 
of threshold 
detection by 
comparison 
of trailing, 
objective and 
leading 
samples. 


Threshold (Diurnal, 
Scintillation and 
Gaseous Attenuation) 


~<— Attenuation, dB 


Control 


where F(@_,, 1, 4,) and F(@_,, 1, @,) are derived 
from system memory (rain model, scaling law and 
simulated system constraints), with the decision 
threshold for @, constrained within +1 V, for 
example. 

If the estimated samples are correct, the tenta- 
tive decision becomes the final decision. If not, d, 
is subtracted from the receive sample b, and the fi- 
nal decision is made on the sign of the difference as 
shown in Figure 9, where the prefiltered bandwidth 
B is specified and the sampling clock rate should 
be high. The clock rate m should be determined by 
trial and error so as to produce a nearly instanta- 
neous check on the measured parameter (rain atten- 
uation versus time). Thus, by setting m large and 
comparing the decision threshold of the objective 
sample with those of the trailing and leading 
samples, we can implement a nearly instantaneous 
tracking of uplink fade. The system memory could 
consist of a random-access memory module where 
the rain model, including scaling law and system 
constraints, is stored and updated. Scintillation, gas- 
eous attenuation and diurnal data are used to set the 
initial threshold. A cost-effective uplink fade-con- 
trol scheme is shown in Figure 10. 
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Power Control Design Pointers 


The following issues should be taken into 
consideration for the uplink power control, specifi- 
cally for future VSAT-type networks: 


¢ The VSAT business is by nature competitive. 
Currently, VSAT terminals range in price from 
$5 ,000-$10,000, depending on the application. 
The time frame for Ka-band VSAT systems is 
perhaps early 21st century, when the unit price is 
expected to be significantly lower. Expensive or 
sophisticated algorithms will not be attractive 
features. 


¢ In most cases, VSAT terminals will be unat- 
tended, so their design should be very reliable 
in operation. 


¢ Present-day VSAT link availability is 0.995, 
which is comparable with terrestrial service 
reliability. 


¢ Solid-state power amplifier capability is limited, 
and hence uplink fade compensation in excess of 
6 to 8 dB may be commercially unattractive in 
most cases. 


¢ Scintillation, gaseous attenuation and diurnal ef- 
fects should be used as a fixed bias term so that 
rain attenuation measurement and control are not 
corrupted by the above secondary effects. It may 
even be better to provide a fixed design margin in 
the uplink EIRP for the secondary effects. 


Conclusion 

The theoretical concept for uplink power con- 
trol discussed here will require further development 
for implementation. Assuming a first-order Markov 
process, a nearly instantaneous tracking of rain fade 
is proposed to reduce prediction error normally as- 
sociated with the process response time lag. 
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Spectrum for Satellite 
Services Allocated 
at WARC-92 


Ann O. Heyward, Communications Systems Branch 
NASA Lewis Research Center 


he 1992 World Administrative Radio 

Conference for Dealing With Frequency 

Allocations in Certain Parts of the Spec- 
trum (WARC~—92), held in Malaga/Torremolinos, 
Spain, concluded on March 3. It produced a num- 
ber of new frequency allocations which will en- 
able new services and systems to be implemented 


Actions taken at WARC-92 — an upgrade to 

the Mobile Satellite Service allocation and an 
allocation to Broadcast Satellite Service (Sound) 
— are the culmination of more than 25 years of 
discussion within the ITU and will permit DBS-R 
services to be implemented worldwide. 


throughout the radio spectrum. Two results of par- 
ticular interest to SATCOM Quarterly readers 
concern the proposal to establish a new General 
Satellite Service (GSS) and consideration of an 
allocation to the Broadcast Satellite Service 
(Sound). 


The United States and other nations proposed 
the definition of a GSS — a new service category 
to accommodate both fixed and mobile satellite 
applications, including Personal Access Satellite 
Service communications. Allocations to this new 
service were proposed at 19.7—20.2/29.5—30 GHz, 
the frequency band to be used by the Advanced 
Communications Technology Satellite to demon- 
strate both fixed and mobile communications capa- 
bilities. Representatives of many nations at the 
WARC were not willing to consider a GSS alloca- 
tion at these frequencies. However, the United 
States achieved an upgrade of the Mobile Satel- 
lite Service allocation to co-primary status with 
the Fixed Satellite Service in the 19.7—20.2/ 
29.5—30 GHz band in ITU Region 2, and in the 
20.1—20.2/29.9-30.0 GHz band in ITU Regions 1 
and 3. The 1992 Conference also recommended 
that an appropriate future WARC consider defining 
a general or multipurpose satellite service and con- 
sider allocating additional spectrum to accommo- 
date growth of multipurpose satellite networks. 
These actions by the Conference represent signifi- 
cant progress in accommodating GSS applications 
in the spectrum. 

Technical efforts in support of an allocation 
to BSS (Sound) in the 500-3000 MHz portion of 
the spectrum have been a pivotal element of the 
joint NASA/Voice of America Direct Broadcast 
Satellite-—Radio (DBS—R) Program. Resolution of 
this issue by the Conference provided an allocation 
to BSS (Sound) at 1452-1492 MHz. Most nations 
will utilize this band for BSS (Sound). However, 
to preserve existing services in this frequency band, 
the United States will use the 2310-2360 MHz 
band; several ITU Region-3 countries and the 
Russian Federation will use the 2535-2655 MHz 
band for this service. These allocations are the cul- 
mination of more than 25 years of discussion of sat- 
ellite sound broadcasting service within the ITU 
and will permit DBS-R services to be implemented 
worldwide. 
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